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j^ , Abstract 

G ' Let n points be chosen randomly and independently in the unit disk. 

'^ ■ "Sylvester's question" concerns the probability pn that they are the 

C , vertices of a convex n-sided polygon. Here we establish the link with 

fi I another problem. We show that for large n this polygon, when suit- 

'""'■ ably parametrized by a function r{(p) of the polar angle <j), satisfies the 

^ . equation of the random acceleration process (RAP), d^r/d</)^ = /(</)), 

K* ' where f is Gaussian noise. On the basis of this relation we derive the 

L^ ■ asymptotic expansion log|?„ = — 2n log n + n log(27r^e^) — con^' ^ + . . ., 

(^ . of which the first two terms agree with a rigorous result due to Barany. 

'nT I The nonanalyticity in n of the third term is a new result. The value 4 

t^^ ' of the exponent follows from recent work on the RAP due to Gyorgyi 

^ '. et al. [Phys. Rev. E 75, 021123 (2007)]. We show that the n-sided 

f^ \ polygon is effectively contained in an annulus of width ~ n~^'^ along 

^ \ the edge of the disk. The distance 5n of closest approach to the edge 

is exponentially distributed with average (2n)~^. 
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Figure 1: An instance of a set of six points chosen randomly in a disk D. The 
convex envelope (dotted) of this set does not pass through the point P and hence 
is a five-sided polygon. This set of points therefore does not contribute to pg{D). 

1 Introduction 



Let n points be chosen independently according to a uniform distribution on 
a disk D C M^. We consider in this work the probability Pn{D) that these 
points are the vertices of a convex n-sided polygon; an example is shown in 
figure [H This question, with the disk D replaced by an arbitrary convex 
domain K C M^, has a long history in mathematics. In 1864 Sylvester [1] 
asked about the value oi p^i^K). 

For a parallelogram S and a triangle T the exact results 
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1 /2n 
n\\n- 



Pn{T) 
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n-l)!3(2n) 
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were shown by Valtr [2l[3]. Since afiine transformations may transform any 
parallelogram into a square and any triangle into an equilateral triangle while 
leaving the distribution uniform, expressions (11. ip do not depend on the 
particular choice of triangle or parallelogram. In the limit of large n, the 
probabilities Pn{K) become very small and we will be interested in their 
asymptotic large-n expansions. Those of Pn{S) and Pn(T) above are given 
by 



logp„(S') = — 2nlogn + nlog(16e^) + C(logn) 
logpn(T) = -2nlogn + nlog{Ye^) + 0{\ogn) 



:i.2) 



in which only the coefficients of the terms linear in n are different. These 
coefficients were subsequently explained by Barany [4|, who derived a general 



result that we will paraphrase as follows. Barany showed that for an arbitrary 
convex domain K of area Ak 



logp„(ir) = -2nlogn + nlog(^ie2^j ^o{n), (1.3) 



in which the meaning of Vk still has to be given. To do so we need the 
concept of the affine length Cc of a convex curve C, defined (see e.g. [1]) as 

Cc= f ^'/'ds, (1.4) 

Jc 

with ds a line element along the curve and k, the local inverse radius of 
curvature. We note that the affine length has the physical dimensionality of a 
(length)^/^. Let now a convex subdomain K' of K have a border dK' of affine 
length Cqk' and let furthermore Cax' attain its maximum for K' = K'^^^. 
Then 

Vk = max Cgx' = Cqk' (1-5) 

is called the affine perimeter of K. This completes the definition of (II. 3p . 

There is no simple way to calculate Vk for general K. When K is a disk 
D of radius Rd, symmetry dictates that -K'max = K = D. Then (11.51) and 
(11.41) with K = Rf^ yield Vk = C,dD = ^ttR^ , so that (11.31) becomes 

hgPniD) = -2n\ogn + n\og{2Tc^e'^) + Tin , (1.6) 

in which the remainder TZn is o{n). This statement about the remainder 
is weaker than that in equation (II. 2p . where the rest term is known to be 
O(logn). It will appear that there is a good reason for this difference. 

In this work we represent the n random points by polar coordinates 
{Rm, ^m) for m = 1, . . . , n in such a way that < $i < $2 <•••<*&« < Stt. 
We define the average radius, i?av, and the scaled deviation from average, 
rm, by 



n 



R,,=n-^J2Rm, i?„ = i?av(l + n-'/V„). (1.7) 

m=l 

For n —>■ 00 with (p = 2'Kmn^^ fixed, r^ becomes a random function r(0) on 
[0,27r]. Our principal result is that in this limit the remainder TZn in (11.61) 
can be cast in the form 

7^„ = log ( exp [ - 2ra^/2 max r(0)])g, (1.8) 

in which (. . .)o denotes the average with respect to all 27r-periodic zero- 
integral solutions r(0) of 



where /(0) is Gaussian noise of autocorrelation 

(/(0)/(0O)o = |[27r5(0-0') -1]. (1.10) 

We arrive at these results by extending an analytic method that was devel- 
oped originally in the context of planar Voronoi tessellations [5|, [6l [7] and 
employed also to study Poisson line tessellations in general and the Crofton 
cell problem in particular [8]. The idea of applying this method here arose 
from the observation that the first two terms in expansion (11.61) are identical 
to those of the expansion of logp^""^, where p^™ is the probability for the 
typical Poisson- Voronoi cell to be n-sided. The work of reference [9] suggests 
that the two problems are related by an inversion of the radial coordinates 
with respect to the unit circle; this is indeed borne out by our analysis. Our 
method consists of a rather intricate but fully exact coordinate transforma- 
tion, followed by an asymptotic expansion which, although nonrigorous, is of 
the kind routinely used in physics to obtain exact results. 

Equation (11.91) is known as the random acceleration process; the function 
r(0) is also referred to as Kolmogorov diffusion or integrated Brownian mo- 
tion. We will briefly review some of the literature on this stochastic process 
in section HI Of particular relevance to us is recent work by Gyorgyi et al. 
(GMOR) [lO]. Upon using their result for the right hand side of (11.81) we 
find that (II. 6p becomes 

\ogPn{D) = -2n logn + nlog(27rV) - 2eo(37r^n)^/^ + . . . , (1.11) 

where eo > is the smallest eigenvalue of a linear eigenvalue problem, the 
only hypothesis being the existence of its solution. An immediate corollary 
is that in the large-n limit the average (II. Sp draws its main contribution from 
polygons that stay within a distance of order ~ u'^^^Rd from the edge of 
the disk D. For n — > cxd the distance of closest approach to the disk is shown 
to be exponentially distributed with average l/{2n). 

The term proportional to n^^^ in (II. lip is a new contribution to the 
answer to Sylvester's question. Since such a nonanalytic term is absent from 
the expansions (II. 2p for the triangle and the square, we are led to ask under 
which conditions such an n^/^ term appears. Barany [1] showed, essentially, 
that in the large-n limit the n points in convex position lie on the curve 
dK'^^^, and that this curve is composed of 

(i) arcs or isolated points that coincide with the domain boundary dK; 

(ii) arcs of parabolas in the interior of K. 
The present study provides strong indication that an n^^^ term occurs when- 
ever dK'^^^ contains at least one arc coinciding with dK, that is, when dK'^^^ 
sticks to the domain boundary over some nonzero angular interval. This is 
obviously the case for the circle, where dK^^^ = dK, but not for the square 
and the triangle, where the limit curve dK'^^^ touches the domain boundary 
dK only in isolated points [^ |3] . 



In section |2] we carry out the exact coordinate transformation. In section 
[3] we perform the large-n expansion and establish relation (11.81) . In section 
14.11 we use GMOR's results to obtain (11.111) . In section 14.21 we show how 
to obtain the power ^ in (11.111) by heuristic arguments. In section 14.31 we 
return to the hypothesis of the existence of the eigenvalue eo; we derive exact 
bounds for the right hand side of (11.81) and show that even without the 
existence hypothesis the conclusion remains valid that 7^„ is nonanalytic in 
n. Section [5] contains the study of the distribution of the distance of closest 
approach of the edge. Section [6] is our conclusion. 

2 Random convex polygon in a disk 

We consider n points Ri, . . . , R„ drawn randomly and independently from 
a uniform distribution on the disk of radius Rn, centered in the origin. We 
ask for the probability Pn{D) that these n points are the vertices of a convex 
polygon. A slightly different probability p* is defined the same way but with 
the additional condition that the polygon enclose the origin. When n gets 
large, the ratio between Pn{D) and p* will tend to unity exponentially rapidly 
with n [H], but p* will be easier to study. We begin by writing its definition. 
In terms of the polar coordinate representation R^ = {Rm, $m) we have 



1 rRo /-Stt 

-^— / R.dRi . . . Rr^dRr, d$i...d$„x(Rl,---,Rn), (2.1) 

<dJ jo Jo 



{ttR 

in which x is the indicator of the subdomain of phase space where the points 
form an origin-enclosing convex ra-sided polygon; the explicit expression of 
X will be discussed in section 12.21 Everywhere below we will scale the 'radii' 
Rm such that Rd = 1. 

We will now subject expression (12.11) to a series of coordinate transforma- 
tions. The final result of these will be equation (j2.22p together with (I2.23p 
and ([22ID. 

2.1 Transformation of variables. I 

In (12. ip we may set one of the angles, say $„, equal to 27r if we compensate 
by an extra factor 27r; and the remaining n — 1 angles may be ordered such 
that 

< $1 < . . . < $„_i < $„ = 27r (2.2) 

if a compensating factor (n — 1)! is introduced. Referring now to figure [2] we 
define the angle differences ^^ between two consecutive vertex vectors Rm-i 
and Km by 

Cm = $m-$m-i, m = 1, . . . , n, (2.3) 



with the convention $o = 0. Then the C,m are positive and satisfy the sum 
rule ^^=1 ^m = 27r. In terms of these variables we can write (12. ip as 

P:= „ ^ / RidR,...RndRr. d6...de„5(^i + ... + en-27r)x. 

TJ" JO Jo 

(2.4) 
Next we define the average radius -Rav and the reduced radii pm by 

1 " 
n ^-^ 

m,=l 

The pm therefore satisfy the sum rule 

1 " 

-Ep- = 1- (2-6) 



n 



m=l 



We now rewrite the integrals on the Rm in (12.41) as 

"1 rl rR^.^ 



[ dRi...dRn = [ dR^^Rl:-' [ ^^ dpi...dpn6(l-n-'\2pr 
Jo Jo Jo ^ m=i 

= / dpi-.-dpn^n-n'^^^Prnj 
•^0 m=l 

X f dR,^Rl:-'lle{R-^'-pm) 

"^0 m=l 



(2.7) 



where 0(^R~^ — Pm) is the Heaviside step function, equal to 1 if R~-^ — pm > 
and to if R~-^ — pm < 0. If (12. 7p is applied to an integrand without i?av 
dependence, then the integration on i?av may be carried out according to 



/■I n 

/ dR..Rl:-'lleiR:;-p 



IIlini<m<n Pm 



'-l-fi-av-'T'av 



2n-l 



f2r2' 



,-1 



max Pm 

l<m<n 



-2n 



(2.8) 



Hence, after we substitute ^Bj in (Q and (I22D in ([23D we get 
{n-l)\ '•2" 



Pn 



riTT' 



-f / dei...de„<5 Ven^-27r 



m=l 



max Pm 

l<m<n 



-2n 



(2.9) 



This completes the conversion of the variables of integration from the R^ 
and $m to the pm and C,m- 
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Figure 2: Figure illustrating the derivation of the convexity condition (j2.1ip . 



2.2 The convexity condition 

We will now find an explicit expression for the convexity condition, enforced 
by the indicator x, in terms of the variables of integration pm and ^m- This 
is most easily done by means of the following geometric argument. In figure 
m the point R is the intersection of the line segments ORm, and Rm-iRm+i- 
The point R^ is in convex position with respect to Rm-i and Rm+i if -R < 
i?m, or equivalently, if 



area(ORm-iRm) + area(ORm+iRm) > area(ORm_iR; 



m+lj 



(2.10) 



The areas are readily expressed in terms of the two angles ^„ 
the three radii Rm-i, Rm, and Rm+i- After division by Rm 
the length scale drops out and ( ]2.10p becomes 



and C,m+i and 

-iRmRm+lRav 



sin^ 

Pm+1 



sin ^^^+1 sin(^^ + ^^+i; 



Pm-l 



Pr, 



m 



1, 



>n. 



(2.11) 



where C,m, and pm are understood to be n-periodic in their index. A similar 
condition, differing only in that all radii are replaced by their inverses, was 
encountered in the case of the Voronoi cells [121 E]; it there expresses the 
condition that the mth point (for m = 1, . . . ,n) contribute a nonzero seg- 
ment to the perimeter of the Voronoi cell. This property is due to a duality 
argument: when C is a convex set containing the origin in its interior, its 



dual (or polar body) C* is the convex set {x G M^; (x, y) < 1 W y E C} 
where (■, ■) is the usual scalar product (see [13]). Formally, when applied to 
a convex n-sided polygon, the transformation provides a new convex n-sided 
polygon such that the projections of the origin onto its edges are given by 
the polar coordinates (-R~^, $m), I < m < n. 

Equations (12. 9p and (12. lip represent an intermediate stage in the trans- 
formation of variables. Equation (12.110 couples the variables of integration 
nontrivially and makes the problem of evaluating (12. 9p a hard one. 

2.3 New angular variables 

In this subsection we introduce the angular variables necessary for the re- 
mainder of our analysis. We refer to figure [31 

Let Sm = {Sm, ^m) bc the projection of the origin onto the line passing 
through Rm-i and R^; the Sm and \l/m may be expressed in terms of the Rm 
and $m- The remaining angles needed in the discussion may then be defined 
in terms of the $m and \l/m- The angle difference r}m between the projection 
vectors S^ and S^+i is 

r/m = ^m+i-^m, m = l,2,...,n, (2.12) 

with the convention \I^„+i = \I^i+27r. The 7]^ satisfy the sum rule ^^^^ ^m = 
27r. Next we define 2n angles between projection and vertex vectors, 

Pm = '^m-^m-l, 'y„, = ^m-^m, 171 = 1, ...,n. (2.13) 

We note that the Pm and 7^ may be negative, as happens when the projection 
Sm falls outside the line segment connecting Rm-i and R^. In any case the 
geometry shows that we must have 

~-<(3m,'ym<-- (2.14) 

For fixed sets of angles ^ = {^m} and rj = {rim} one may still jointly rotate 
the vertex vectors R^ with respect to the projection vectors S^, as this 
modifies only the relative angles Pm and 7^ (see figure [3]) between the two 
sets. We may select any one of these relative angles and call it 'the' angle of 
rotation, since it will determine all others. We will select Pi as this special 
degree of freedom and express the remaining Pm and 7^ as 

m—l 

Pm = Pi-^i^i-Ve), m = 2, . . . , n, 
£=1 

m— 1 

7m = -Pi + ^i^i-Ve)+^m, m = 1, . . . , n, (2.15) 

1=1 

8 




m—1 



Figure 3: Figure illustrating the definition of the angles ^m-, Vm, f^m-, and 7^- 

Now observe that the reduced radii pm satisfy the relations 

^.£2£|i, ™.l....,„, (2.16) 

Pm COS (3m 

where po = pn- This ratio (12.161) is the same as the one encountered in the 
Voronoi and Crofton cell problems [6l[8], except for an interchange of Pm o^nd 
7m- It follows that one may express the pm exclusively in terms of the angles 
by solving them from fl2.16p together with fl2.6p . 

We return now to the convexity condition (12.111) . We may use (12.161) 
to eliminate pm±i from (12. lip in favor of pm, which subsequently divides 
out. After some trigonometry (12. lip appears to reduce to the condition 
tan7m + tan/5m+i > 0. Because of (I2.14p this is equivalent to 'jm + Pm+i > 0. 
Still using (12.150 we see that 

r]m>0, m = l,2,...,n, (2.17) 

is, finally, the condition for the n points to be in convex position. 



2.4 Periodicity in the polar angle 

Let us next define the function G by 

n 

2nG(i,mM TT cos 7m 



n 



\ cos/3„ 



f2.18) 



9 



with the /3m and 7^ given by fl2.15p . This function appears when we relate 
Pn to po by taking the product of fl2.16p on m = 1, . . . ,n. Because of the 
periodicity condition po = pn the angle of rotation Pi cannot be arbitrary: it 
must have the special value Pi = P<t{^,ri) that is the solution of 



Gi^,v;P*) 



0. 



(2.19) 



This equation is invariant under the interchange of Pm and 7^ and is there- 
fore identical to the one that appeared in the Voronoi and Crofton cell prob- 
lems. Indeed, as has already been said, the dual body of a convex n-sided 
polygon is a n-sided polygon such that the sets of angles {^^ ; 1 < m < n} 
and {rjm \ 1 < m < n} (as well as {Pm ; 1 < m < n} and {7™ ; 1 < -m < n}) 
have been exchanged. It was shown in reference [H] that the solution /3*(^, if) 
of f l2.19p . subject to fl2.14p . exists and is unique if and only if the ^m and r]m 
satisfy the condition 



m— 1 



m—1 



l<m<n 



max ^(6 - Tjt) + im - min V(^^ - rj^ 

<m.<r> ' ^ \ <m.<n ' ^ 



< TT. 



(2.20) 



We will let 6(^,77) denote the indicator function, equal to zero or to unity, 
of the domain in phase space where f l2.20p is fulfilled. 



2.5 Transformation of variables. II 

Finally we carry out the transformation that eliminates the pm occurring in 
(12. 9p in favor of the rjm defined in (I2.12p . The details of this transformation 
step have been described in [6] and the appendix of [Hj, and we will not 
reproduce them here. The result is an integral on the two sets of variables ^ 
and rj. It is described conveniently with the aid of a 'zeroth order' probability 
distribution V^^\^,ri) defined as 



V'^'K^.r,) 



1 



m=l 



n n 

<^(5Z^m-27r)5(5^r7„-27r), (2.21) 



m=l 



m=l 



where TV = (27r)'^"^^/[(2ri — l)!(n — 1)!] is the appropriate normalization 
constant. The average of any quantity X with respect to V^'^^ will be denoted 
as {X)o. The transformation then leads to 



p:=pi°^(0e-%, 



pi°) 



(Stt 



2\n 



47r2 {2n)\ 
in which G is defined following (I2.20p and V is given by 



(2.22) 



1 



G'i^,v;P*) 



n 



Pm sin ^„ 



^^ ^m cos Pm cos -f„ 



max Pm 

. l<m<n 



-2n 



(2.23) 
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where the prime on G denotes differentiation with respect to its last argu- 
ment. The factor 1/^m in the product in fl2.23p is compensated by the ^m 
in the product in (12.211) . This has been so arranged in order that exp(— V) 
do not diverge for any of the ^^ tending to zero. We note in passing that 
Pn is the same zeroth order probabihty that occurs in the Voronoi problem 
[51 [13]. Equations (I2.2ip - fl2.23p are fully exact and are at the basis of all 
further developments. 

3 Large-n limit and Random Acceleration Pro- 
cess 

The initial problem (12. ip has been transformed into the evaluation of the 
average (Be~^)o in (12.221) . This, however, is still a formidable problem. 
/^From here on we will be able to proceed only by making a large-n expansion. 
This is the subject of this section. 

The average involves only angles and we make the hypothesis, to be con- 
firmed by the calculation, that for large n they all scale with negative powers 
of n. To begin with, we study the scaling that follows from V^'^' {$, , r]) . At a 
later stage we will then discuss how this scaling is modified by the presence 
of the integrand Be"^ in f lT^ . 

3.1 Scaling in the large-n limit. I 

We start with a preliminary. We note that according to V^^^ the ^m are 
independent apart from the sum rule constraint represented by the delta 
function; and a similar statement holds for the rim- Furthermore, the ^m 
and rji are mutually fully independent. More precisely, let {X^, m > 1} and 
{Yi; / > 1} be two independent sequences of i.i.d. random variables such that 
Yi is exponentially distributed with mean 1 and X^ is distributed with law 
u{x) = 4xe~^^, X > 0. For a fixed n, the variable ^^ (resp. rji), I < m,l <n, 
is equal in law to i2TcXm)/{Y.l=i^k) (resp. to C^irYi) / {Y.l=iYk)) ■ Indeed, 
let us average any positive bounded test function h{rii, . . . ,r]n) on [0,27r]"' 
with respect to the set {Yj; Z > 1}. A direct change of variables provides 

= /27r)"-i" / dr/i...dr7n/i(//i,...,r7„)e I ^r/fc-27r J , (3.1) 

which shows that the average with respect to the exponentially distributed 
i.i.d. Yfc is the same as the average weighted with V^^\ The same method 



11 



applies to the Xm and ^m- We now have the following consequence. Let 
us consider the limit n ^ oo with k a fixed integer and (mi, m2, . . . , mk) a 
not necessarily fixed subset of {1, 2, . . . , n}. Then in this limit the marginal 
probability distribution of the vector (?2/27r)(?7m^, . . . ,77^,^) converges to the 
probability distribution of (l^j, . . . , l^t.), i-c, to a product of A; exponentials. 
An analogous statement holds for the ^^ and X^. 

Let us define the scaled zero- average variables 5xm and 5ym by 

n~'^5xm = Cm- 27™"\ 

n~^6ym = r]m-27Tn'^, m = l,. . . ,n. (3.2) 

They satisfy the sum rules J2m=i ^^"^ ~ Ylm=i ^Vm = and one readily 
calculates their variances and covariances, 

{5Xi5Xm)o = —-(Sim ), 

{SyeSy^)o = — -—[Sim ), 

n + 1 V n/ 

{Sxt5ym)o = 0, i,m = l,...,n, (3.3) 

which will be needed later. The necessary calculations may be carried out 
directly via (12.211) or with the use of the explicit realization of C,m and rim as 
functions of Xi, ■ ■ ■ , X„ and Yi, ■ ■ ■ , y„. Equations (13. 3p show that 6xm and 
Sym are of order n° in the large n limit. 
We define Um by 

W^Um = |(/?m-7m) 
m—1 

= (3i - ^{Ci - Vi) - l^m , m = l,...,n, (3.4) 



where to obtain the second line we used (I2.15p . It follows from (13.41) together 
with the scaling of the 6xm and 6ym found in 13.21 that in the large-n limit 
the (3m and 7^, being sums of independent random variables, have Gaussian 
distributions of average vr/n and of width ~n~^/^. The Um therefore remain 
of order n° when n gets large [15] . We define r^ by 



Pm = l+n ^Tm, m = l,...,n. (3.5) 

/^From the scaling of jSm and jm together with (I2.16P and (13. 5p it follows that 
the Tm are of order n". Due to the sum rule on the pm in (12. 6p they satisfy 
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3.2 Second order recursion 

We are now ready to perform the large-n expansion of equation fl2.16p . 
Pm/Pm-i = cos 7m/ cos /5m, using the scahng of the preceding subsection. 
This yields a recursion relation which in terms of the scaled variables r^, 
Um, and 6xm takes the form 

rm-rm-i + ... = |(/3m-7m) + --- (3-6) 

where the dots represent terms of higher order. Using that Pm + Im = ^m 
and Prn — lm = 2n^^^'^Um we get from this in the large-n limit 

n n 

where now on both sides we have kept only the lowest order. Here and 
throughout the remainder it is understood that all quantities are n-periodic 
in m. In (13. 7p the first and the second term on the right hand side come from 
the average and the random part of the angle ^mi respectively. Both are of 
order n~^, but there is a difference. To see this we imagine to sum fl3.7p on an 
m interval of order n^ for some 5 > 0. Since Um has long-range correlations 
we may consider it as effectively constant over this interval. Then, because 
of the 5xmi the second term will become a sum of independent zero-average 
random variables and hence its contribution will be of relative order n~^/^ 
with respect to that coming from the first term. We may therefore neglect 
the SxmUm term in (13.70 when our purpose is to study the variation of r^ on 
a scale that increases as a power of n. Taking now second order differences 
we find the recursion 

3 

rm_i -2rm + rm+i = 27rn"2/^, (3.8) 

with the right hand member defined by 

fm = -\{5Xm-'25ym + 5Xm+l)- (3.9) 

This quantity satisfies the sum rule ^^=1 fm = 0. 

3.3 Second order differential equation 

We may pass to the following continuum description. We define the polar 
angle = 27cm/n and the functions r(0) = r^ and dr/d0 = u{(f)) = Um- 
In the large-n limit r becomes a continuous function of a continuous vari- 
able. For the second order difference on the left hand side of (13. 8p we get 
Tm-i — 2rm + fm+i -^ (27r/n)^d^r/d0^. On the right hand side we put 
fm -^ 27111^ 2 f {(f)) . This converts (13.80 into the second order differential 
equation 

^^ = /(0), O<0<2vr, (3.10) 
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valid for angle differences |0 — 0'| on the scale n"^"*"*^, that is, large on the 
scale of the discrete index m. In the same way as in the previous subsection 
we may imagine to sum fl3.10p on an m interval of length ~ n"^. In the large-n 
limit ^^l^Q fm is a Gaussian variable centered at zero, whose correlation 

with similar variables on neighboring intervals is small. Hence on scales ^ n^ 
the right hand side of (13.101) becomes Gaussian noise. Its autocorrelation 
function follows from definition fl3.9p and the correlations (13.31) . namely 

(/(0)/(0O)o = |[27r5(0-0') -1]. (3.11) 

The noise is white except for the extra —1 term inside the brackets. This 
term leaves all Fourier components of / with nonzero wavenumber unaffected: 
they behave as under white Gaussian noise. Only the zero-wavenumber com- 
ponent of / is exceptional. By integrating (13.111) on and 0' one finds that 
Jq '^d0 /(0) = 0, in agreement with the sum rule in the last line of section [3^ 
Now the differential equation (13.101) shows that instead of having this sum 
rule on / one may equivalently impose on r that its derivative be periodic, 
r'(0) = r'(27r). The functions r(0) that we must deal with are therefore, in 
summary, those that satisfy equation (13.101) with boundary conditions 

r(0) = r(27r), r'(0) = r'(27r) (3.12) 



and that, due to the relation in the last line of section 13.11 obey moreover 
the zero-integral constraint j„^d(j)r{(j)) = 0. For this class of functions the 
noise / is effectively white. 

In the statistical physics literature equation (I3.10p represents what is 
called the Random Acceleration Process (RAP). Some time ago this process 
appeared [6j (but without being named explicitly) in a similar way in the 
study of many-sided Voronoi cells. In section l^?T] of this work we will see that, 
through the connection that we have discovered here, known facts about the 
RAP provide further answers to Sylvester's question. 



n 



3.4 Large-n expansion of j9* 

We return to the evaluation of the probability p* given by (I2.22|) . 

P: = pi°^(0e-%, Pi°^ = ^^, (3-13) 

in which V is given by (I2.23p . All operations that led to these equations 
were exact for any finite n. We will now see what happens if we perform a 
large-n expansion assuming the scaling of section 13.11 In the large-n limit 
the condition imposed by 0(^, rj) in (12.221) will be satisfied with a probability 
that goes exponentially fast to unity; hence we have 9 ~ 1 and 

P:^pf(e-^)o, (3.14) 
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where the ~ sign denotes an equahty up to corrections that are exponentially 
small in n. We consider now the n ^ oo limit of V. In (12.231) we have 
G'{^,rj; 13^) = 1 + 0{n~^) [H]. Expanding the other factors for small angles 
we get 

-V 

e = exp 



1 v^ , / \2 



+ ... 



n 

m=l 



exp 



— I "d0 {r\<P) - i—f^ - 2ni/2 /0) 

27r io ^ M0W 0<<^<27r ^^^ 

+ ( max r(0)) +...], (3.15) 

where the dots indicate terms of higher order in the angles. As n — >■ oo with 
the scaling discussed in section [3TT1 the dot terms in the exponential in (13.151) 
tend to zero; the first term (which is a sum of n terms compensated by a 
factor n~^) and the third term stay of order vP] and the second term is of 
order n^/^. 

3.5 Scaling in the large-n limit. II 

As shown by (I3.14p . expression (I3.15P has to be submitted to an average, 
denoted (. . .)o, with respect to all solutions r(0) of the random acceleration 
process satisfying the constraints stated in section 13.31 The law for r(0) 
follows, via equation (13.101) . from the law for /(0), which in turn follows, via 
(13.91) . from the law (12.211) for the ^rn and r^i. When n is large, the second term 
in the exponential in (13.151) will suppress large excursions of r(0) from zero 
in a severe but as yet quantitatively unknown way. Let us therefore suppose 
that the average (e~'^)o draws its principal contribution from those r(0) that 
stay within a tube of a width Wn ~ n'"" with an as yet unknown a. Hence 
the effectively contributing r(0) will scale as 

r(0) ~ n~°, max r(0) ~ ra"" (3.16) 

for some as yet unknown a > 0. We must ask ourselves, first of all, if we 
are allowed to have the scaling already introduced above followed by this 
new scaling. In fact we are, because it only restricts us more narrowly to 
the center of the realm where the first scaling holds. The same results are 
arrived at if the full combined scaling is used from the start. The reason for 
the procedure we chose is that this is the best way of showing the connection 
with the Random Acceleration Process, of which finally only a subprocess 
plays a role, namely the one consisting of trajectories that stay within a tube. 
Secondly, if indeed the main contributing r(0) scale with a negative power 
of n, then the first and the third term in (I3.15P are negligible with respect 
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to the second one in a large-n expansion. So finally the problem becomes to 
calculate the average (e~'^)o with 

(e~^)o~ (expr-2n^/2 max rU)])^, (3.17) 

where in the last line the average is with respect to the process r(0) defined 
by fl3.10p and (13.111) . In equation (13.171) the only reference to n is the one 
explicitly visible in the exponential on the right hand side. Hence here n 
has become an "external" parameter coupling to the maximum value of an n 
independent random acceleration process. Finally, the remainder 7^„ in the 
expansion of the convexity probability Pn{D) in a disk is related to (I3.17P by 

7^„ = log(e-^)o. (3.18) 

The connection, in the large-ra limit, between Sylvester's question and the 
Random Acceleration Process is the principal achievement of this paper. We 
will now see that, when combined with what is known about this process, it 
leads to the results announced in the introduction. 



4 Properties of the Random Acceleration Pro- 
cess 

4.1 The work by Gyorgyi et al. 

The random acceleration process is the k = 2 member of the wider class of 
equations 

d^f/d^'^ = /(0) (4.1) 

with fc = 1, 2, . . ., and where /(</>) is white Gaussian noise. Because of their 
relevance in physics and in applied statistics, many examples of such pro- 
cesses have been studied. 

In early work on the k = 1 case, where f(0) is Brownian motion, Foltin 
et al. [16] analyzed the distribution of the square width of r{4>). Burkhardt 
[T7] considered the k = 2 problem in a half-space. More recently, Majumdar 
and Comtet [HI [19], were interested in the distribution of the "maximum 
height" Tmax, defined as the maximum of f(0) relative to its interval average, 

1 r'i-K 



rmax = max f (0) - — / d0f(0). (4.2) 

O<0<27r Zn Jq 

Using powerful path integral methods, Majumdar and Comtet showed that 
for k = 1 and periodic boundary conditions the probability law of the max- 
imum height is the Airy distribution. They pointed out the importance of 
this problem as an instance of extreme value statistics of a set of strongly 
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correlated random variables. The equation with general k was studied very 
recently by Gyorgyi et al. [10] (GMOR), who also focused on the distribution 
-P(^max) of the maximum height (14.21) . The authors of reference [10] calculate 
this distribution for the class of functions f(0) that satisfy (for k = 2) the 
periodicity conditions 

f(0)=f(27r), r{0)=f{2n). (4.3) 

Hence, after subtraction of the interval average, this is the same class that 
appears in the averages (. . .)o in sections 13.41 and 13. 5[ This establishes the 
applicability of GMOR's results to the problem of section [31 

By extending Majumdar and Comtet's [TS] methods GMOR found an 
expression for the distribution P(rinax) of rmax- Still more recently, Burkhardt 
et al. considered the distribution of the maximum relative to the initial value 

m- 

Expression (I3.17P is the Laplace transform with Laplace variable 2n^/^ 
of the maximum height distribution of the random function r{(j)). No math- 
ematically rigorous results for this quantity seem to exist; however, closely 
related properties, involving in particular the maximum of the absolute value 
\r{(f))\, were studied, e.g., by Khoshnevisan and Shi |21j and more recently 
by Chen and Li [22]. The Laplace transformed maximum height distribution 
did figure among the many properties of the RAP studied by GMOR jTD] . 
These authors express this quantity by means of a trace formula involving 
the second order differential operator 

19^ d 

which is the generator of the Markov process (r(0), ^(0)). They consider the 
linear eigenvalue problem 

Cip{r,u) = —eilj{r,u), < r < oo, — cso < u < oo, 

V'(0,m) = 0, u>0, (4.5) 

and assume that it has a well-defined solution with lowest eigenvalue eo. 
When converted to our notation [23], the result of GMOR, obtained essen- 
tially by a scaling of all variables in the trace formula, is that in the limit of 
large n 

- log(e-^)o - 2eo(37r^n)^, n ^ oo. (4.6) 

With (14.61) we have obtained our main result, equation (11.111) of section [H 

4.2 Heuristic argument for the power 1/5 

We will derive the power | found by GMOR by a heuristic scaling argument. 
This is not meant as additional evidence for this result, but as a physical way 
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Wn~n 




Figure 4: Heavy solid arc: edge of the disk of radius Rd = 1- Solid curve: the 
convex polygon R{<j)) in the limit n — > oo. Dashed arcs: annular tube of width 
Wn ~ n~^'^. The two dotted transverse line segments mark the ends of a typical 
tube section of length A^ between two consecutive 'reflections' of the polygon at 
the inner and the outer tube walls; the section length scales as A^ ~ n~^'^. The 
distance of closest approach of the edge scales as 6n ~ n~^. 



of understanding it. We focus on the n dependence of (e~^)o. If we agree 
to consider the expression for V in the exponential in ( J3.17P as an 'energy', 
then what we will present is basically an energy- versus-entropy argument. 

We begin by supposing that due to the factor n^/^ included in the 'energy' 
in (13.171) the contributions to (e~'^)o will come essentially from trajectories 
r(0) with a maximum less than some small but as yet unknown width Wn- 
These contributing trajectories will typically be confined to a tube of width 

2Wn. 

We refer now to figure HI We imagine the tube divided into sections such 
that inside each section the trajectory evolves freely, but at the sections' ends 
it is 'reflected' by the tube walls. Let A^ be the angular interval of a typical 
section (and hence A„i?av is its length). We may estimate the n dependence 
of A„ and Wn as follows. If u{(j)) = dr/d0 is the 'radial speed', then (13.101) 
may be rewritten du/d(f) = f{4>). In an angular interval A„ this speed will 
change by 6u = J, " d(j)'f{(f)') ~ An , so that the typical radial speed inside 

the tube will be of order \6u\ ~ A™ . It follows that in this interval the 



XnSu 



A 



3/2 



This means 



,3/2 



value r(0) of the radius itself will change by Sr 

that we have typically \r\ ~ maxo<(/,<27r ?"(0) ^ Wn ^ Kl^ ■ Hence a typical 
contributing trajectory r(0) will contribute to the average on the right hand 
side of (I3.17P an amount exp(— est x n^/"^ x An ). This is our result for the 
'energy' factor. 

We now estimate which fraction of all trajectories is 'contributing' in the 
above sense, that is, stays within the tube of width 2wn- This confinement 
costs a constant amount of entropy per section, due to the 'reflection' at its 
ends, the trajectory being essentially free everywhere else. Since the number 
of sections is 27r/An, our result for the 'entropy' factor is exp(— est x 2-n\~^). 



Upon multiplying the energy and entropy estimates we arrive at the 
heuristic estimate for (e~'^)o. Taking logarithms yields 

- log(e-^)o ~ C,\-' + C2nV2A3/2, ^ ^ oo, (4.7) 

where Ci, C2 > 0. In (14.71) the n dependence of the section length A„ is still 
unspecified. By setting A„ ~ n^° and minimizing the right hand side with 
respect to a we find a = |, whence 

A„ ~ n-l/^ Wn ~ n-3Ao. (4.8) 

We may return now to the original unsealed variables. By transposing (11.71) 
to continuum notation we see that the radius of the polygon, R{(f)) = Rav[l + 
n~^/^r(0)], runs through a tube of width Wn = Rs^n'^^'^Wn ~ n~^^^, as 
depicted in figure HI 

4.3 Exact bounds 

The mathematical methods used in section[2]are fully exact. Those of section 
[3] are basically a systematic expansion in negative powers of n; although we 
did not rigorously prove its validity, it has been obtained by methods that 
commonly lead to exact results. The argument of section l^?Tl due to GMOR, 
involves an assumption of a different kind, namely the existence of a lowest- 
energy solution of (14. 4p and (14. 5|) with an eigenvalue eo. We have no reasons 
to doubt that this solution exists. However, we wish to point out here that 
without this existence assumption one may derive by elementary methods a 
weaker but still very meaningful statement. We have relegated the proof to 
Appendix A. We there show by straightforward methods starting from (13.171) 
that we have the bounds 



1 . . , \T. . 1 



Cgrae < — log(e )o < 04^4, n —^ 00, (4.9) 

with positive constants C4 and cq. 

Hence, with or without the existence assumption of eo, it is established in 
either case that the large-n expansion of log(e~'^)o depends nonanalytically 
on n. Consequently, because of (13.181) and (II. 6p . the expansion of logpn{D) 
for a disk contains this same nonanalytic term and is clearly distinct from 
the expansions (11.21) for squares and triangles. 



5 Closest approach of the edge 

Let us call 6n the distance of closest approach between the perimeter of the 
polygon and the edge of the unit disk, that is, 

6n = I — max Rm- (5.1) 

l<m<ra 
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Given the results of the preceding sections, it is now fairly easy to study this 
quantity. 

We now first turn to the study of the distribution of -Rav, which we will 
call P(i?av)- We have 

P(i?av) = ^%^, (5.2) 

Pn 

where the denominator p* is given by integral (12. ip and the numerator p(-Rav) 
by this same integral except that we put primes on the variables of integration 
and insert a factor 5{R'^^ — -Rav)- Evaluation of the numerator proceeds as it 
did for the denominator in section 12.11 Instead of (12.81) we have to calculate 

/ di?l,i?li"-i5(i?l,-i?av) = i?L"-^e([ max p^y -i?av). 

J n. l<rn<n 

(5.3) 
It then follows that 

P(^-) = 7' .v\ ' (5-4) 

(Be ^)o 

in which exp(— V^,^^) is the same expression as (I2.23P except that the factor 

[maxi<m<n Pm]"^" has been replaced with 2ni?^""^ ([maxi<„<„ p^]"^ - -Rav) • 

Expression (15. 4p for P(-Rav) is still fully exact. 

We cannot go beyond this point unless we make again a large-n expansion. 

We define a variable (iav by 

i?av = 1 - ri"'/'4v, (5.5) 

where we tacitly suppose that (iav is at most of order r/' but possibly smaller. 
Using in (15. ID that Rm = RsArPm, employing (15. 5p . and expanding, we obtain 
between 6n and (iav the relationship 

6n = n"-^/^ [dav - max Tm] + . . . , (5.6) 

'- l<m<n -' 

where the dots indicate terms of higher order. 

Using expression (13.51) for p^ and expanding as in section 13.41 we get an 
expression for P(-Rav) in terms of the variable (iav, 

P(i?av) = 2n(e~V(exp[— 5^(r^-«^)-2nV24^-rfL + ...l) • (5.7) 



n 

m=l 







For the same reasons as in section 13.41 we neglect now the sum and the 
d1^ term in the exponential in (15. 7p . Finally, using (15.60 . we eliminate the 
variable (iav from (15. 4p in favor of 5„ whose distribution we will call p{6n). We 
get in the limit n — > cx) the final result of this section, namely the probability 
distribution p of the distance of closest approach of the edge, 

p(5„) = 2ne-2"'5", 5n>0. (5.8) 

Not only does this show that 6n scales as 1/n, but also that it is exponentially 
distributed. 
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6 Conclusion 

In this work we considered Sylvester's question: given n points chosen ran- 
domly from a uniform distribution on a disk, what is the probability that 
they are the vertices of a convex n-sided polygon? For n large this probability 
becomes very small and one may consider the asymptotic expansion of logp„. 
The first two terms of the expansion were known; they are proportional to 
nlogn and to n. In this work we establish, first of all, a relation between 
Sylvester's question and the Random Acceleration Process. We then show 
that the third term in the expansion (in absolute value) is asymptotically 
bounded from above and below by ~ n^/^ and ~ n^/^, respectively, so that 
it must be nonanalytic in n. If one accepts the hypothesis underlying the 
work of Gyorgyi et ai, which we easily do, it follows that this term must be 
proportional to ~ n^^^ with a well-defined coefficient given in (11. lip . Along 
with this expansion of p„ we harvest a variety of results concerning the most 
probable way that the n points are distributed along the edge of the disk. 

The subject has not been fully exhausted. Remaining questions in the 
present work concern, for example, correlations between the angles. Besides, 
a natural generalization would be to consider more general convex domains 
K and in particular convex polygons. In this context, the work [24J provides 
precise results on convex chains for deriving exact distributional results on 
the number of vertices of the convex envelope. Extensions of this work could 
deal with the crossover phenomenon that must exist when the convex do- 
main K varies between finite-sided polygons and the disk: for example, what 
happens when one considers n points randomly chosen in the interior of a 
regular M„-gon when n and M„ tend to infinity in a specified way? Another 
extension would be the study of the convex envelope of iV„ randomly chosen 
points knowing that this envelope is n-sided, again in the limit of n and Nn 
tending to infinity. However, we leave these and other matters for future 
investigation. 



A Bounds for (e ^)o as n ^ oc 

In order to derive upper and lower bounds on the average (e~'^)o we begin by 
expressing this quantity in terms of the Fourier transforms of the variables 
involved. 



A.l Fourier transforms 

We define the Fourier transform 

-. n n 

h = r E e'"'^'"/"/™ , r, = n"^ V e^^'^^/V^ , (A.l) 
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where, if for convenience we take n odd, g = 0, ±1, ±2, . . . , ±(|n — \)- The 
sum rules imply that fa = tq = 0. In Fourier language recursion (13.81) 

becomes 

7r2 . 1 . 

^'? = - ^2si^2^ /g - --ih^ 5^0. (A.2) 

n " 

where the ~ sign indicates the limit n —^ oo aX fixed q. This amounts to 
neglecting higher orders in g/n, which are small if in agreement with our 
preceding discussion we restrict ourselves to g on a scale ~ n^"*^, that is, 
to spatial distances which in units of the index m scale as ~ n^ . lYiom. 
definition (jA.ip together with (13.91) we have furthermore 

(/J,')o = |V.',o, g,gVo. (A.3) 

The fq are distributed according to a marginal distribution P^^\f) of P*^°^(^, r]). 
In the limit n — i> oo all fq become Gaussian distributed and therefore P*-''-' is 
given by 



,„ V 37r 

This distribution may be used for averaging quantities that depend essentially 
on the long wavelength properties of the process. 

It will be easier to deal with real, as opposed to complex, quantities. 
Choosing a convenient normalization we define the radial and angular com- 
ponents Fq and ojq of f±q by 

f±q= \^Fq{cOSUJq± ism UJq), g > 0. (A. 5) 

Letting as before = 27r?7i/n, we have that in terms of these Fourier trans- 
forms the process r(0) becomes 

oo 

r(0) = ^ e-^'?'^ fq = -VQ Y1 -i ^°^ (^<^ - ^'?) ■ (^-6) 

q q=l ^ 



Substituting flA.6P in (13.171) and transforming to the Fq and Uq as new vari- 



ables of integration we get 

"Ztt oo , oo 






?2 
_ _ Q 

q=l '"^ q=l " q=l 



X exp 



oo 



— 2v6n max >^ ^ cos (g0 — cjg) . (A. 7) 



Q<6<2tt ^ — ^ a 



In the right hand side of (]A.7p the only n dependence that comes in is the 
one explicitly exhibited in the argument of the exponential. We note that the 
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Fourier expression (]A.6|) shows that r(0), due to the factor 1/q^ in the sum 
on g, is precisely a quantity of the type that receives its main contribution 
from small q values and whose properties may therefore be calculated by 
averaging with respect to (lA.4p . The feature that makes this problem hard 
to solve exactly is the maximum that is required in flA.7p . In what follows 
we will obtain upper and lower bounds for (e~'^)o. 

A. 2 Two inequalities 

We recall that the argument of the second exponential in flA.7p is equal 
to — 2Y^maxo<0<27r ''"(0)- Upper and lower bounds for (e~'^)o are based on 
bounds for this maximum expressed in terms of the Fourier amplitudes of 
r(0). These bounds, valid asymptotically for n large, read 

2 

— < max rid)) < M, (A. 8) 

where 

g=l ^ g=l ^ 

the Fq are defined by flA.6l) . and in which it is assumed that both sums on 
q converge. The upper bound given by (]A.8P - (]A.9P is obvious and leads to 
a lower bound for (e~''^)o derived in Appendix IA.3I The lower bound given 
by (lA.Sp - IJA.Op is proved in Appendix lA. 4. It it leads to an upper bound for 
(e~^)o which is derived in Appendix IA.4.2[ 

A. 3 Lower bound for (e ^)o 

In this section of the appendix we abbreviate 



N = 2VQn. (A.IO) 

We will find a lower bound for (e~''^)o valid in the limit of large A^. We replace 
the maximum in ( JA.7P by its upper bound given in ( lA.Sp . Since this bound 
is independent of the Ug, the integrals on these variables reduce to a factor 
unity. The integrals on the Fg factorize and we get 

(e"^)o > n / 2a;da;exp (-x'^-Nq~^xj = exp -^J^iqN'i) (A.ll) 

g=l -^0 q=l 

where 

poo 

J^[k) = -log / 2xdx exp(-x^ - k^'^x). (A.12) 

Jo 
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This function is positive for all k > and it is integrable since it behaves as 
JF(k) ~ 41ogK~^ for K ^ and as J^(/t) ~ ly^/t"^ pour k — > oo. The sum 

on q in ( JA.lip makes the argument of JF vary by steps of spacing A^~2 and 
we may therefore in the limit N ^ oo replace it by an integral, which gives 



e"^)o > exp -A^2 / dKj^(K 
^ Jo 

= exp [ — C4n4] , 



(A.13) 



where in the last step we used the relation between A^ and n given in (JA.IOI) 
and where C4 is a positive constant. 



A. 4 Upper bound for (e ^)o 

In order to prove the upper bound for (e~'^)o, we first need to prove the lower 
bound on r(0) stated in (1A.8I) . 



A. 4.1 Lower bound for m.a,Xo<^<2n r{(j)) 

We prove the following property. 

Property. Let f{x) be a function on [0, 27r] such that 

"277 



/ dxfix) = 0, 



dxf{x) = 27^a^ 



min fix) = — M, 



where a and M are given positive constants. Then 

4^2 \ 1/2 



max f{x) > \M 



0<x<27r 

In the limit M ^ a this yields 



M2 



a 



max fix) > -— 

0<a<27r ^ ^ ~ M 



2 
\APJ 



i + o(^) 



Proof. We set 



/(x) = /+(x)-/_(x) 



(A. 14) 

(A.15) 
(A.16) 



(A.17) 



(A.18) 
(A.19) 
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in which f±{x) 


= maxjO, ±f{x)\. In terms of these Eqs. flA.14)-(A.16D be- 


come 










/■27r /■27r 

/ dxf+{x) - / dx/_(x) = 
Jo Jo 


= 0, 


(A.20) 




/ dxflix)+ / dxflix) -- 
Jo Jo 


= 2rta^ 


(A.21) 




max f (x) = 


= M. 


(A.22) 


We abbreviate 










m = max f\(x). 

0<x<2n-'^^ ' 




(A.23) 



We now have the estimate 



1 f^"" 
m > — / dxf+(x) 

1 f'^'^ 
= 2^1 ^^^-(^^' 



(A.24) 



where in the second step we used (1A.20|) . Next 



27r 



dxf^(x) > — T / dxf'^(x) 



2tt 



^[i^"' - [" <ixfl(x)] 



^ |[--'"^]' 



(A.25) 



where in the first step we used (1A.22I) . in the second step ( 1A.21I) . and in the 
third step flX23l) . 

We now combine flA.24p and flA.25P to get 



m > 



cr^ — m^ 



M 



(A.26) 



By solving the associated second-order equation for m and using the def- 
inition flA.231) we are led directly to flA.17p . This completes the proof of 
property flA.lTp . 

Application. We apply inequality flA.17p to the function r{(f)) given in 
( ]A.6|) . Its M satisfies M < Mq with Mq given in equation (\AM and its a^ is 
given in that same equation. Since (as is implied by the arguments of section 
IA.4.2P for n — i> oo the ratio ct/Mq tends to zero, we may replace (1A.17P by 
(lA. ISp and obtain the desired lower bound stated in (lA.Sp . 
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A. 4. 2 Proof of upper bound for (e ^)o 

We will next find an upper bound for (e~''^)o valid in the limit of large n. 
We replace the maximum in (1A.7P by its lower bound given in (lA.Sp . The 
integrals on the variables uOq reduce again to a factor unity. The integrals on 
the Xq do not factorize in this case. In order to make them do so we introduce 
integral representations and write 

(e-^)o < 






n2F,dF,exp[-^F|-liv5:-j/$: 

q=l q=l q=l ^ q=l 

Jo J~oo ^TT Jo 1 ^ \=1 ^ ^ q=l ^/^ q=l 



d/^ r 

J —oc 




N 'log / 2xdxi-x^ J 

^U Jo ^ ^m' 


d/xexp 


\-^'0] 


fj-- 


\ r( ^M 



2 I'S 



in which 



^1(2;) = 2^1ogag(2;), al{z) = l + — 



9=1 

00 



6^2(5,2;) = -y;iog/ 2|/di/exp(-y2 + _!f^^y (A.28) 

We will need the large z behavior of Qi{z). Upon scaling q = kz* in ( JA.28P 
we find 



/•oo 

^1(2;) ^ ao-24, ao= dfi;log(l + K"^), 

Jo 



z ^ 00. 



(A.29) 



In order to study Q2{,s,z) we expand the integrand in its definition (1A.28P 
in a power series in is, do the y integration term by term, and expand the 
logarithm in powers of is. This gives 



Q2{S,Z) 



g=l k=0 



k\ 



q^Uq 



V^Ai{z)is - 1(1 + TT)A2{z){isf + ... 



where 



Adz) 



00 _. 



A;=l,2, 



(A.30) 
(A31) 



(A.27) 



26 



After again scaling q = kz^, we get by the same method as above 

Ak{z)=akz'^--2''[1 +...], ak= d/€(l + K^)-^^ z -^ oo, (A.32) 

Jo 

where the dots stand for a series in powers of 2;" 2 . Let us denote by J„ (yU, g-) 
the integral on s in the last line of (lA.27p . Substitution of the above expansion 
of Q2 yields 



"'"(/^' 2u) 



^ ds 
-— exp 



27r 






-h'^(l + 7i){a2 



+ 



+ 0{{g] %^ 



(A.33) 



in which the dots come from expansion (1A.32I) . Let us scale s = (^) *t. We 
then get 



^"(-^) - &[^ 



X exp 



^^^^-^)'-"^"^(^)'^'-'(' + ")"^'' 



27r 



'1 + 7r)a2 






/^ 



2/1 



1 fN\h^ 



2(l + 7r)a2 



+ ... 
(A.34) 



where the dots stand for terms proportional to negative powers of N/{2fi). 
It is now possible to see how we should scale /i with A^. We will introduce 
the scaled variable z/ defined by 2/i = z/^A^~3, where the fourth power on u is 
just a matter of convenience. Substituting flA.34p in flA.271) . using the large 
z expansion flA.29p . and neglecting terms that vanish as A^ ^ cxd we then get 



(e-^)o < est X A^3 / dz/e 



-N^g(u) 



with 



9iu) = c,u-' + ^' 



1 _l\2 

2 — aiTTii/ 2 j 



(A.35) 



(A.36) 



2(l + 7r)a2 

The function g{u) increases as ~ u^^ for z/ ^ and as ~ z^ for z/ ^ cxo, 
and hence has a minimum value for some z/ = z^min- Upon doing the integral 
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(IA.35|) by steepest descent we get 

(e~^)o < est X A^2 exp [ - Nig{v^ 

= est X n4 exp [ — Cgne] , (A. 37) 

where Cg is a positive constant. 

Acknowledgments 

HJH kindly thanks Zoltan Racz for a useful exchange of correspondence. 

References 

[1] J.J. Sylvester, Problem 1491, The Educational Times (April 1864), Lon- 
don. 

[2] P. Valtr, Discrete Comput. Geom. 13, 637 (1995). 

[3] P. Valtr, Combinatorica 16, 567 (1996). 

[4] I. Barany, The Annals of Probability 27, 2020 (1999). 

[5] H.J. Hilhorst, J. Stat. Mech. L02003 (2005). 

[6] H.J. Hilhorst, J. Stat. Mech. P09005 (2005). 

[7] H.J. Hilhorst, Eur. Phys. J. B: Proceedings of Statphys 23, Genoa, Italy, 
July 9-13, 2007. Available online (2008). 

[8] H.J. Hilhorst and P. Calka, to appear in J. Stat. Phys. (2008); 
larXiv:0 802. 18691 

[9] P. Calka and T. Schreiber, Ann. Probab. 33, 1625 (2005). 

[10] G. Gyorgyi, N.R. Moloney, K. Ozogany, and Z. Racz, Phys. Rev. E 75, 
021123 (2007). 

[11] A polygon not enclosing the origin is confined to a half-disk. The asser- 
tion then follows from a result due to J. G. Wendel {Math. Scand. 11, 
109 (1962)). It can also be seen as a consequence of applying 01. 3p to a 
half-disk. 

[12] P. Calka, Adv. m Appl. Probab. 35, 863 (2003). 

[13] R. Schneider, Gonvex bodies: the Brunn- Minkowski theory (Cambridge 
University Press, Cambridge, 1993) 

28 



[14 
[15 

[16 

[17; 

[18 
[19 
[20 

[21 
[22 
[23 



H.J. Hilhorst, J. Phys. A 40, 2615 (2007). 

Sums such as im = I3m + 1m and 77^ = 7^ + Pm+i require caution, since 
the leading ~n~^/^ behavior of the two terms cancels and the result is 
of order n~^. 

G. Foltin, K. Oerding, Z. Racz, R. Workman, and R.K.P Zia, 
Phys. Rev. E 50, R639 (1994). 

T.W. Burkhardt, J. Phys. A 26, L1157 (1993). 

S. Majumdar and A. Comtet, Phys. Rev. Lett. 92, 225501 (2004). 

S. Majumdar and A. Comtet, J. Stat. Phys. 119, 777 (2005). 

T.W. Burkhardt, G. Gyorgyi, N.R. Moloney, and Z. Racz, Phys. Rev. E 
76, 041119 (2007). 

D. Khoshnevisan and Z. Shi, Trans. Amer. Math. Soc. 350, 4253 (1998). 

X. Chen and W.V. Li, The Annals of ProhahiUty 31, 1052 (2003). 



Our r(0) is a factor v37r larger than the function h studied by GMOR 



[24] C. Buchta, Mathematika 53, 247 (2006). 



29 



